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Leptonic/hadronic models for electromagnetic emission in 
microquasars: the case of GX 339-4 
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ABSTRACT 

We present a general self-consistent lepto/hadronic jet model for the non-thermal 
electromagnetic emission of microquasars. The model is applied to the low-mass mi- 
croquasar (LMMQ) GX 339-4 and predicts its high-energy features. We assume that 
both leptons and hadrons are accelerated up to relativistic energies by diffusive shock 
acceleration, and calculate their contribution to the electromagnetic spectrum through 
all main radiative processes. The radiative contribution of secondary particles (pions, 
muons and electron-positron pairs) is included. We use a set of simultaneous obser- 
vations in radio and X-rays to constrain the model parameters and find the best fit 
to the data. We obtain different spectral energy distributions that can explain the 
observations, and make predictions for the high-energy emission. Observations with 
gamma-ray instruments like Fermi can be used to test the model and determine the 
proton content of the jets. Finally, we estimate the positron injection in the surround- 
ing medium. Our findings support the suggested association between LMMQs and the 
observed distribution of the 511 keV line flux observed by INTEGRAL. 

Key words: gamma-rays: theory - radiation mechanisms: non-thermal - X-rays: 
binaries - X-rays: individual: GX 339-4. 



1 INTRODUCTION 

The low-mass microquasar GX 339-4 was discovered in 1972 
by the satellite OSO-7 (Markert, Canizares & Clark 1973). 
Since then, it has been extensively observed at all wave- 
lengths from radio to X-rays, and detected in all the ca- 
nonical spectral states of X-ray binaries. Little is known, 
however, about the remaining characteristics of the binary 
system. 

Based on modulations in the optical photometry, 
Callanan et al. (1992) inferred an orbital period of 14.8 hs, 
later confirmed by Buxton & Vennes (2003) using optical 
spectroscopy. Further optical spectroscopic measurements 
and analysis of long-term X-ray light curves showed no ev- 
idence of this modulation, revealing instead a periodicity 
of ~ 1.75 days (Hynes et al. 2003; Levine & Corbet 2006). 
The first estimates of the distance to GX 339-4 placed the 
system at d ~ 1.3 — 4 kpc, see Zdziarski et al. (1998) and 
references therein. This result was later revised by Zdziarski 
et al. (2004), who concluded that the minimum distance lay 
in the range 6.7 kpc < d m i n < 9.4 kpc. They favoured a lo- 
cation in the galactic bulge at ~ 8 kpc. However, a study of 
absorption lines performed by Hynes et al. (2004) suggests 
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that a location in the far side of the Galaxy at a distance 
d > 15 kpc cannot be completely ruled out. 

The emission in the optical band is dominated by the 
accretion flow (Imamura et al. 1990), preventing direct ob- 
servation of the secondary star even when the system is go- 
ing through the very low X-ray luminosity state. The first 
detection of the star was made by Hynes et al. (2003) during 
the X-ray outburst observed in 2002. The mass and spectral 
type of the donor star have not been firmly established yet. 
According to Hynes et al. (2004), an orbital period of ~1.7 
days implies a low density for the companion of ~ 0.06 g 
cm~ 3 . This probably corresponds to a low- mass subgiant of 
spectral type G or F, depending on the assumed distance. 
Following this idea, Munoz-D arias, Casares & Martmez-Pais 
(2008) suggested that the star is a "stripped-giant" , in which 
mass loss is due to the burning of an hydrogen shell. In 
this model, the mass of the secondary must be in the range 
0.166M© < M 2 < 1.1M©. From this result, the authors con- 
strained the mass of the compact object to be Mbh > 6Mq 
or even Mbh > 8.6M0, for a mass of the secondary near 
the lower or upper limit, respectively. These values strongly 
support the idea that the compact object is a black hole (see 
also Hynes et al. 2003). 

GX 339-4 has been observed in radio, infrared, opti- 
cal and X-ray wavelengths, sometimes simultaneously or 
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quasi-simultaneously (Hannikainen et al. 1998; Wilms et 
al. 1999; Nowak, Wilms & Dove 2002; Homan et al. 2005). 
The source goes through all the spectral states of X-ray 
binaries: low-hard, high-soft, very high state, intermediate 
state and quiescence (McClintock & Remillard 2006). It fre- 
quently displays outbursts associated with state transitions, 
episodes during which the X-ray luminosity can reach peaks 
of Lx = 10 37-38 erg s" 1 for an assumed distance of 6 kpc 
(Homan et al. 2005; Yu et al. 2007). It was after the X-ray 
outburst of 2002 that Gallo et al. (2004) imaged for the first 
time a relativistic radio jet on ~ 10 3 AU scales in the sys- 
tem (see also Corbel et al. 2000). The detection of the jet 
confirmed that GX 339-4 is a microquasar. 

The hard X-ray emission during the low-hard state of 
microquasars and black hole binaries is generally thought 
to have its origin in a hot corona that surrounds the black 
hole. The observed power-law spectrum is explained through 
Compton up-scattering of accretion disc photons by hot elec- 
trons in the corona. Some of these up-scattered photons in 
turn excite iron nuclei in the disc material, giving rise to 
the appearance of a Ka line at ~ 6.5 keV, superimposed on 
the power-law continuum; for evidence supporting the pres- 
ence of the Fe Ka line in GX 339-4, see for example Dunn 
et al. (2008). However, Corbel et al. (2003) showed that in 
GX 339-4 the radio and X-ray fluxes are tightly correlated, 
Jradio oc F x ' 7 ■ This suggests that the emission in both bands 
might have a common origin in synchrotron radiation pro- 
duced by non-thermal electrons in the jet, and not in the 
corona (Corbel & Fender 2002; Corbel et al. 2003). This idea 
was explored by Markoff et al. (2003) and Markoff, Nowak 
& Wilms (2005), who applied a jet model to fit the observa- 
tions. They showed that synchrotron radiation of relativistic 
electrons in the base of a jet can explain both the radio and 
X-ray spectra and reproduce the observed correlation. This 
model, however, is purely leptonic: the contribution to the 
radiation output of relativistic protons in the jet is not taken 
into account. 

In this work we present a lepto-hadronic model for the 
broadband electromagnetic spectrum of GX 339-4. In our 
model, protons and electrons are accelerated up to relativis- 
tic energies. The hadronic radiative contribution extends 
well into the gamma-ray band, since typical maximum pro- 
ton energies are much higher than those of electrons. We 
expect that our predictions can be tested in the near fu- 
ture with the data collected by instruments like the Fermi 
gamma-ray satellite and by atmospheric Cherenkov arrays 
like the High Energy Stereoscopic System II (HESS II). 

Along with electromagnetic radiation, the creation of 
electron-positron pairs is a necessary result of relativis- 
tic particle interactions. They are injected, for example, 
through photon-photon annihilation and as a by-product of 
hadronic interactions. 

Recently, measurements carried out with the INTE- 
GRAL Spectrometer (SPI) instrument of the INTEGRAL 
satellite have allowed to complete a detailed map of an ex- 
tended region of emission line at 511 keV in the Galaxy 
(Weidenspointner et al. 2008). These observations confirm 
the diffuse (rather than point-like) distribution of the line, 
with bright emission around the galactic center (flux ~ 10~ 3 
cm~ 2 s" 1 ), and a clear, asymmetric (toward negative longi- 
tude values) disc component. The disc total flux is ~ 1/5 of 



the bulge flux. The disc emission is detected up to scales of 
~ 20° from the galactic center. 

Different types of positron sources have been suggested 
in the literature, including pulsars, the massive black hole 
at the galactic center, microquasars, nucleosynthesis events, 
and extended processes like cosmic ray nuclear reactions and 
dark matter decay. The fact that both bulge and disc emis- 
sion are clear, seems to disfavour the latter two possibilities 
as the main positron sources. Microquasars, in particular 
those with a low-mass donor star, seem to be a particularly 
appealing possibility, given the spatial correlation of the line 
emission with the overall distribution of low-mass X-ray bi- 
naries in the Galaxy (Weidenspointner et al. 2008) and en- 
ergetic arguments (Guessoum, Jean, & Prantzos 2006). 

Up to the present, none of the microquasar models avail- 
able presents quantitative predictions of the positron pro- 
duction rate, although in Romero & Vila (2008, 2009) the 
radiative output of positrons was taken into account. The re- 
cent findings described above makes it timely an exploration 
of the possibilities for electron-positron production in self- 
consistent models for microquasar jets. In the present paper 
we have devoted a section to estimate the positron produc- 
tion rate in systems like GX 339-4. 

The article is organized as follows: in Section f2] we 
present the basis of the model, describing in detail the hy- 
pothesis made about the geometry and energetics of the sys- 
tem. We also describe the calculation of the relativistic par- 
ticle distributions (for primary and secondary particles) and 
the radiative processes. In Section [3] we discuss the obser- 
vational data sets and the constraints they impose on the 
model parameters. We also present the best-fitting spectral 
energy distributions obtained and analyse the results. Our 
estimations for the positron production rate are described 
in Section |4] Finally, in Section [5] we summarise our results 
and we compare them with those of previous works. We close 
discussing the predictions for the very high-energy emission. 



2 MODEL 

2.1 Basic picture 

The jet model applied here is based on the model developed 
in Romero & Vila (2008). A sketch of the system is shown 
in Figure [T] The jet is supposed to be conical, and launched 
at a distance zo = 50i? g from the compact object, where 
R g = GA/bh/c 2 is the gravitational radius of the black hole. 
The initial radius of the jet is taken to be a fraction \ °f the 
value of the injection distance, ro = X z o- Following Falcke & 
Biermann (1995), we relate the jet power to the Eddington 
luminosity of the black hole as 

Lj c t = — (ftet^Eddj (1) 

where L E dd ~ 1.3 x 1O 38 (A/ BH /M ) erg s" 1 . The factor 
1/2 accounts for the existence of a counterjet of equal power. 

At an arbitrary distance z from the compact object, the 
total energy budget of the jet can be roughly divided into 
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Figure 1. Scheme of the jet. It is launched at a distance zo of the 
black hole (BH) , and acceleration of particles takes place further 
away at z acc . The width of the acceleration region is indicated as 
Az. The jet axis makes an angle (j> with the line of sight. 



magnetic energy, bulk kinetic energy and particle internal 
energjQ 



Lj et (z) w L B (z) + L k i n {z) + L int (z). 



(2) 



If the plasma is ejected by some kind of magnetocentrifugal 
mechanism, we can assume that at the base of the jet the 
entire jet power is in the form of magnetic energy, Lj ct ~ 
Lb(zo). This allows to estimate the value of the magnetic 
field Bo = B(zo) by equipartition between the magnetic 
energy density and the total energy density, 



8n 



-'jet 



(3) 



where Vj c t is the bulk velocity of the outflow. After the 
launching point the magnetic field decreases with the dis- 
tance z to the compact object as 



»(.)-»(?)■ 



with 1 ^ m 2 (Krolik 1999). Magnetic energy is then 
converted into bulk kinetic energy (and a fraction of this is 
dissipated at shocks and converted into internal energy of 
relativistic particles, see below), in such a way that equipar- 
tition holds only at the base of the jet. 

A few remarks are in order here. In the region z < z aC c 
the outflow is a Poynting-dominated flux, formed by a ther- 
mal plasma that accelerates as it gets energy from the or- 
dered magnetic field. Thermal electrons and ions in the "pre- 
acceleration" region move approximately with the jet bulk 
velocity and have an energy ~ Tjctmc 2 . The plasma is neu- 
tral, but almost all the bulk kinetic energy of the jet is car- 
ried by the ions. No significant electron synchrotron radia- 
tion is expected from this region. At z ~ z acc shocks perturb 
the field, and a fraction of the energy associated with the 
bulk motion of the plasma is converted into internal en- 
ergy of relativistic non-thermal particles. In this region syn- 
chrotron cooling becomes efficient, especially for electrons. 



1 Strictly, Lj ct > -^jet(^) 
ation. 



since 



part of Lj e t is dissipated as radi- 



This situation is different from the case of high-mass mi- 
croquasars, where "cold" relativistic electrons can cool by 
inverse Compton interactions with the stellar photon field. 
A similar scenario arises in binary systems formed by a pul- 
sar and a luminous Be star, such as PSR B1259-63/SS2883, 
see Khangulyan et al. (2007). 

Most of the bulk kinetic energy of the jet is carried by a 
thermal plasma. If shock waves propagate through some re- 
gion of the outflow, the suprathermal tail of the Maxwellian 
particle distribution can be accelerated up to relativistic en- 
ergies by diffusion across the shock front (see, for instance, 
Drury 1983 and references therein). The physical conditions 
for an efficient acceleration are not clear. For the plasma 
to be mechanically compressible and allow the formation of 
shocks, the magnetic energy density Ub = B 2 /8tt must be 
in sub-equipartition with the bulk kinetic energy density [7k 
of the plasma (see Komissarov et al. 2007 for a discussion 
on this topic). Therefore, the acceleration region has to be 
located at a distance z aC c from the black hole such that 

Ub(Zzcc) < C4(Zacc)- (5) 

The kinetic energy density of the jet can be written as 



Uk = n(z)E p " 



(6) 



where n(z) is the cold particle density, and Ep ln is the rel- 
ativistic kinetic energy of a proton that moves with the jet 
bulk velocity, 



Ep ln = (Tjet - l)m p c. 



(7) 



The presence of shocks may not suffice to accelerate 
particles efficiently. According to Gaisser (1990), for diffu- 
sive shock acceleration to work, the ram pressure in the ac- 
celeration region must dominate over the magnetic pressure. 
This condition can be written as 



(8) 



where Um is the internal matter energy density. For a cold- 
proton-dominated jet, Um can be calculated as in Bosch- 
Ramon, Romero & Paredes (2006), 



U m = 



7rrf et ujetm P 



(9) 



(10) 



Here Mj e t is the jet mass flow, 

and -Ep ln is the classical kinetic energy of a thermal proton, 

(11) 



£,kin 1 2 

E p = 2 m P v p- 

The mean velocity of the particles was taken to be equal 
to the lateral expansion velocity of the jet, v p — w e x P = X v iet , 
that is of the order of the speed of sound in the plasma. Con- 
dition J8J is stronger than ([5} , in the sense that if the former 
is fulfilled, so is the latter. In either case, the location of the 
innermost acceleration region can be determined demanding 
that the appropriate condition is satisfied, 



Ub — pU(k, m ), 
with p < 1. 



(12) 
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We assume that the jet is dynamically dominated by 
cold matter. The total power L re i injected in relativistic par- 
ticles is only a small fraction of the jet power, 



(13) 



This energy is shared between relativistic protons and lep- 
tons, L rc i = L p + L e . We relate the energy budget of both 
species as 



(14) 



The parameter a remains free in our model but, as we are 
interested in jets with a relevant hadronic content, we kept 
a ^ 1 throughout. 

2.2 Fundamental equations 

2.2.1 Particle distributions 

The differential spectrum of particles accelerated by diffu- 
sion through shock waves is a power-law in energy, and the 
particle injection function can be written as 



Q(E) = QoE~ 



(15) 



The normalization constant Qo is calculated from the total 
injected power, 



d 3 r 



EQ (StP) (E) dE, 



(16) 



where V is the volume of the acceleration region. Particles 
can gain energy up to a certain value iJ max for which the 
sum of the cooling rates 

(17) 



1 



E dt 



over all the relevant processes of energy loss i, equals the 
acceleration rate (Aharonian 2004), 

ecB 



t: 



(18) 



The parameter r\ < 1 characterises the efficiency of the ac- 
celeration mechanism. 

The steady-state energy distributions of relativistic par- 
ticles N (E) is then calculated solving the transport equa- 
tion in the one-zone approximation in the jet co-moving 
reference frame (where the distributions are supposed to be 
isotropic) , 



d_ 

dE 



dE 



dt 



N(E) 



N(E) 
T 



= Q(E). 



(19) 



The time T is the characteristic time of "catastrophic" non- 
radiative particle losses: processes through which particles 
are removed from the system, for example by escape or de- 
cay. As protons and electrons are stable particles, T is simply 
equal to the escape time-scale from the acceleration region, 



Az 

Wjet 



(20) 



Equation (|19[) neglects effects of convection and diffu- 
sion, and therefore it is only valid in a thin region of the jet. 
Its exact analytical solution can be found, for example, in 
Khangulyan et al. (2007). 



2.2.2 Radiative processes 

We consider several mechanisms of interaction of the rela- 
tivistic particles with the magnetic field, photons and matter 
in the jet. Both protons and leptons interact with the mag- 
netic field emitting synchrotron radiation. Electrons also ra- 
diate by relativistic Bremsstrahlung as they are accelerated 
in the mean electrostatic field of the ions in the jet plasma, 
and by inverse Compton scattering (IC) against the syn- 
chrotron photon field (Synchrotron Self Compton process, 
SSC). Interaction of protons with the synchrotron field (py) 
proceeds through photopair production 
„+ 



p + 7 — > p + e + e 
and photomeson production 
p + 7 — > p + an + b (n + + 7r~) , 
p + 7 — > n + 7r + + air + b (ty + + ty 



(21) 

(22) 
(23) 



An approximation for the cross section for photopair pro- 
duction is given in Maximon (1968). If e — x'm e c 2 is the 
photon energy in the proton rest frame, then for x' < 2 



1.2 x 



+ 



mb, 



(24) 



whereas for x' > 2 

<r e = 0.58x [3.1 In 2a/ -8.07 + ...] mb. (25) 

In the case of photomeson interactions, the cross section av 
can be roughly approximated as a step function (Atoyan & 
Dermer 2003), 

f 340 ^ibarn 200 MeV e ^ 500 MeV 
ov(e') ?s <^ 

[ 120^ibarn e > 500 MeV, 

Relativistic protons also produce pions by interaction with 
the jet matter field through proton-proton (pp) collisions, 



p + p 

p + p 
and 
p + p 



• p + p + an + b (tt + + n ) , 

• p + n + 7r + + air + b (ty + + 7r~) 



n + n + 



2?r + + an° + b (tt h 



(26) 
(27) 

(28) 



For high proton energies, the probability of creation of the 
three pion species is almost the same. The total inelasticity 
of the process is ~ 0.5; most of the kinetic energy lost by the 
proton is carried away by only one or two leading pions. The 
inelastic pp cross section can be accurately approximated as 
(Kelner, Aharonian & Bugayov 2006) 



Cinel (Ep) 



34.3 



l. 88m JpL_ +0 .251n 2 Ep 



ITeV 



ITeV 



1 - 



Eth 
E„ 



mb, (29) 



where E t ^ — 1.22 GeV is the proton threshold energy for ir° 
production. The pion multiplicities a and b have power-law 
dependences on the relativistic proton energy, a,b oc Ep K 
with n ~ 1/4 (Mannheim & Schlickeiser 1994). 
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The gamma-ray output from proton-proton and photomeson 
interactions is due to the decay of neutral pions, 



2 7 . 



(30) 



Besides radiative cooling, particles also lose energy 
through adiabatic losses, since they exert work on the jet 
boundary surface. The expression for the adiabatic cooling 
rate and for all the radiative mechanisms of energy loss are 
compiled in Romero & Vila (2008) and Reynoso & Romero 
(2009). 

We calculate the contribution of each radiative pro- 
cess to the luminosity. For synchrotron radiation, relativistic 
Bremsstrahlung, and inverse Compton scattering (in both 
Thomson and Klein-Nishina regimes) we use the formu- 
lae given in Blumenthal & Gould (1970). The gamma-ray 
luminosities from proton-proton and photomeson interac- 
tions are calculated following Kelner, Aharonian & Bugayov 
(2006) and Kelner & Aharonian (2008), respectively. 

All the luminosities, except those from process of in- 
teractions with matter, are better calculated in the jet co- 
moving reference frame. The corresponding luminosity in 
the observer frame is then obtained applying the Doppler 
boosting factor D((j>), 



D= [r jct (l -# ct cos0)] 1 



(31) 



where <f) is the viewing angle (the angle between the jet axis 
and the line of sight) and /3j e t = Wjet/c. Denoting with primes 
the magnitudes in the co-moving frame, for an approaching 
jet the luminosity transforms as 

L 7 (£ 7 ) = D 2 L' 1 (E! y ), (32) 
whereas the photon energy in the observer frame is 

E-y = DE!y. (33) 



2.2.3 Secondary particles 

Charged pions, muons and electron-positron pairs produced 
in proton-proton and proton-photon interactions can also 
contribute to the electromagnetic emission spectrum. If the 
magnetic field and/or the photon and matter fields in the jet 
are strong enough, and fi can significantly cool through 
synchrotron radiation, inverse Compton scattering and even 
pion-proton collisions (np) before decaying. Charged pions 
decay creating muons and neutrinos, 

-> n~ (34) 

7r+^M + + iV (35) 

Muons in turn decay yielding electrons, positrons and 
more neutrinos, 



e + Vu + v e 



e + Vn + v e 



(36) 
(37) 



The pion injection function Q 7T ± for pions created in pp 
collisions is given in Kelner, Aharonian & Bugayov (2006); 
for pions injected through photomeson interactions we ap- 
ply the 5-functional approximation, see Atoyan & Dermer 
(2003). From Q 7T ±, we calculate the muon injection function 
Qu± using the formalism presented in Lipari, Lusignoli & 
Meloni (2007). 



Electron-positron pairs created by muon decay and 
those directly injected through photopair interactions are 
also a source of electromagnetic radiation. We calculate the 
pair injection function from muon decays following Schlickei- 
ser (2002). To estimate the pair output Q e ± from photopair 
production, we apply the formulae given in Chodorowski, 
Zdziarski & Sikora (1992) and Miicke et al. (2000). Finally, 
there is another source of pair injection due to photon- 
photon (77) annihilation, 



7 + 7 ■ 



(38) 



The corresponding pair source function is given in Bottcher 
& Schlickeiser (1997). The process of 77 annihilation is also 
a photon sink, and can eventually modify the production 
spectrum. Its effect on the escape of photons is discussed in 
the next section. 

The steady state distribution of secondary pions, muons 
and electron-positron pairs is calculated solving the trans- 
port equation (|19|) . in the same way as for the primary par- 



ticles. In the case of n ± and fj, , the decay time must be 
included in T along with the escape time, 



f- 1 4-t' 1 

CSC ~ L c\f}r.' 



(39) 



In the jet co-moving reference frame, the decay time 
t dcc is given by 



idee — T(7T,n) 7(jt,m)' 



(40) 



(7r,M)/ m (ir.M) c2 ^ S Lorentz factor of the par- 



Here 7 (lrjM) = E { 
tides and T(^ m j is the decay time in their rest frame, 
TV = 2.6 x 10"' 8 s and r u = 2.2 x 10" 6 s. 



2.2.4 Attenuation of the production spectrum 

As remarked above, gamma-ray photons can annihilate 
against low-energy radiation to create electron-positron 
pairs. For a gamma ray of energy E 1 to interact with a 
photon of energy e, their energies must exceed the threshold 



e£ 7 (1 



(41) 



The angle 8 is that between the momenta of the collid- 
ing photons. According to equation (I41|l . TeV gamma rays 
can be strongly absorbed by infrared photons produced, for 
example, by electron synchrotron radiation in the jet. Thus, 
this process of internal absorption can modify the shape 
of the high-energy region of the production spectrum. The 
probability of absorption of a gamma ray can be quantified 
through the opacity t 77 , defined as 



Tjj(Ej) — ^ 



cr 77 (_E 7 , e, u) n(e) x 

(l-u)dudedL (42) 



Here cr 77 is the angle-averaged cross section (Gould & 
Schreder 1967), n(e) is the distribution of target photons, 
u — cos 9, and / is the length of the path traversed by the 
gamma ray from the emission site to the observer. The shape 
of the modified luminosity L 7 can be estimated from the pri- 
mary production spectrum as 



Lj {E~f) — exp (~t 77 ) Lj (E~f) . 



(43) 
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This provides a zeroth-order approximation to the final 
emission spectrum, mainly indicating those energy ranges 
where the emission is suppressed. A more refined treatment 
of the problem requires solving a set of coupled equations 
to calculate self-consistently the equilibrium distribution of 
particles and photons, or to solve the equations of an electro- 
magnetic cascade, if it develops. However, the strong mag- 
netic field suppresses IC cascades through fast cooling by the 
synchrotron channel (see Khangulyan, Aharonian & Bosch- 
Ramon 2008). Cascades are expected to play a more im- 
portant role in high-mass X-ray binaries, see for example 
Orellana et al. (2007). 



3 GENERAL RESULTS 
3.1 Application to GX 339-4 

GX 339-4 was extensively observed simultaneously in ra- 
dio and X-rays during the low-hard state (LHS) in 1997, 
1999 and 2002. For some of these observations, simultaneous 
near infrared (NIR) and optical data are also available. The 
1997 and 1999 radio observations were carried out with the 
Australia Telescope Compact Array (ATCA) and the Mo- 
longlo Observatory Synthesis Telescope (MOST), and are 
described in detail in Corbel et al. (2000). The radio, NIR 
and optical data from 2002 is presented in Homan et al. 
(2005) . The X-ray data were collected with the Rossi X-ray 
Timing Explorer (RXTE), and are compiled in Wilms et al. 
(1999), Nowak, Wilms & Dove (2002), Corbel et al. (2003) 
and Homan et al. (2005). We refer the reader to these works 
for model assumptions and other details of the data extrac- 
tion in each case. No further data reductions were performed 
for this paper Additional information on the particular 
data sets used here is presented in Table [T] 

We placed the injection point of the jet at a distance 
zo — 50i?g ~ 4.5 x 10 7 cm from the black hole, and fixed the 
ratio x = r o/zo = 0.1 (Romero & Vila 2008). This gives a jet 
half-opening angle of m 6° . The position of the acceleration 
region was determined from equation ([8]), demanding sub- 
equipartition of the magnetic energy density with respect to 
the jet internal matter energy density at z — z acc . A second 
scenario in which z acc is fixed from equation ((5J was also 
considered. We constrained the ratio of energy densities to 
be p < 1, but otherwise it was left as a free parameter. In any 
case the acceleration/emission region was taken to be thin, 
Az — 4z acc , as required for the one-zone approximation to 
be valid. We adopted several decay prescriptions for B(z), 
with m =1.2, 1.5, 1.8 and 2. 

Obsl, Obs2 and Obs5 correspond to the end of the low- 
hard state, when the source was highly luminous. Assuming 
a conservative value of d = 6 kpc for the distance, the ob- 
served X-ray fluxes yield luminosities of up to Lx ~ 10 37 

2 Calibration and data reduction algorithms have been updated 
since the data presented here were reduced. Reprocessing the data 
might result in changes in the slope of the X-ray spectrum. How- 
ever, we do not intend to perform detailed fits to the spectrum 
but to show that the observations can be accounted for by a lepto- 
hadronic jet model, as an alternative to purely leptonic models. 
Using the same data allows to compare the results of our model 
with those of previous works. 



erg s _1 . This places some constraints on the value of the pa- 
rameters that determine the energetics in our model. Only a 
small fraction of the jet power is carried by relativistic par- 
ticles, otherwise the outflow could not be confined; we fixed 
<7roi = 0.1 in equation (113[) . In a model with equipartition 
between hadrons and leptons (a — 1), half of this energy is 
given to relativistic electrons. If the observed X-ray flux is 
due to electron synchrotron radiation, this implies at least 
a total jet power Lj ct ~ 2 x 10 38 erg s _1 . This is a signifi- 
cant fraction of the Eddington luminosity of a black hole of 
Mbh = &M Q , L E dd « 7.8 x 10 38 erg s" 1 . If part of the ac- 
cretion power is radiated outside the jet and part advected 
onto the black hole, the accretion rate required to account 
for the observations must be very near the Eddington limit. 

An accretion model that could apply to powerful sources 
or high luminosity states has been proposed by Bogovalov & 
Kelner (2005). They showed that, along with the standard 
thin disc solution of Shakura & Sunyaev (1973), there exists 
another accretion regime in which the disc is radiatively very 
inefficient, even for high accretion rates. In this solution, 
known as the "dissipationless disc model", a magnetized 
plasma falls onto a central object. The plasma is attached 
to the magnetic field, and angular momentum is removed 
from the system not by viscosity effects, but it is carried 
away by matter itself (see also the Advection-Dominated 
Inflow-Outflow Solution model of Blandford & Begelman 
1999, for early ideas regarding this accretion regime). In 
fact, the model predicts that the mass advection rate van- 
ishes at r = 0, and all the infalling matter is ejected. In this 
way, most of the accretion power could be directly channeled 
into the jets. This model could account for the observations 
of very powerful jets and low-luminosity discs in extreme 
systems such as SS433 or M87. Other radiatively ineffi- 
cient models include Advection-Dominated Accretion Flows 
(ADAFs, Narayan & Yi 1995) and Magnetically-Dominated 
Accretion Flows (MDAFs, Fragile & Meier 2009). These 
models, however, are more suitable for low accretion rates. 

GX 339-4 was also detected during the low luminosity 
phase of the low-hard state in 1999 (Obs3 and Obs4). The 
typical X-ray luminosity is Lx ~ 10 34 erg s _1 . Applying 
the same energetic considerations as above, the minimum 
jet power required is now Lj e t ~ 2 x 10 35 erg s _1 , a fraction 
gj ct « 3 x 10~ 4 of the Eddington luminosity of the black 
hole. 

The observed spectrum in the X-ray band is quite hard, 
Lx oc E~ p with p « 0.3. If the X-rays originate in electron 
synchrotron radiation, from the observed slope it is possible 
to estimate the spectral index 5 of the steady-state parent 
particle distribution, N oc E~ s . They are related as 




This yields S « 2.4. Since particles cool, the index S is not 
the same as that of the electron injection function, Q oc E~ a . 
In particular 8 = a + 1 in the case of dominant synchrotron 
losses. The particle injection spectrum must therefore be 
quite hard, with a power- law index smaller than the typically 
assumed a — 2.0—2.2 predicted by the theory of acceleration 
in strong, non-relativistic shocks. Here we decided to fix a = 
1.5, consistent with relativistic shock acceleration (Stecker, 
Baring & Summer ling 2007). 

The values of the relevant parameters of the model are 
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Table 1. Observational data used in the model fits. 



Observation 


Date 




X-ray flux 




Radio flux 






(y/m/d) 


(10"' 


3 erg cm -2 s - 1 ) 




(mJy) 








3-9 keV 


9-20 keV 


20-200 keV 0.8 GHz 


4.8 GHz 


8.6 GHz 


Obsl 


1997/02/03 


1.06 


1.02 


4.95 7.0 




9.1 


Obs2 


1999/04/02 


0.49 


0.48 


2.75 


4.8 


5.1 


Obs3 


1999/06/25 


0.059 


0.052 


<0.29 


0.14 


0.34 


Obs4 


1999/08/28( a > 


0.037 


<0.01 


<0.17 




0.35 






X-ray flux 


IR/optical magnitudes 


Radio flux 








(10 - 9 erg cm -2 s 






(mJy) 








3 - 300 keV 


m H 


mj my 


4.8 GHz 




Obs5 


2002/03/22 


15.1 


11.7 


14.1 15.6 


13.3 





M Radio data from 1999/09/01. 

X-ray data were collected with the Proportional Counter Array (PCA) and High Energy X-ray Timing 
Experiment (HXTE) instruments of the RXTE satellite. The HXTE measurements are normalized to PCA 
flux levels. Optical and IR photometry was obtained with the Yale-AURA-Lisbon-Ohio State (YALO) 
telescope. The radio flux density at 0.8 GHz was obtained with MOST, and at 8.6 GHz and 4.8 GHz with 
ATCA. From Corbel et al.(2000), Nowak, Wilms & Dove (2002), Corbel et al. (2003) and Homan et al. 
(2005). 



Table 2. Values of the relevant parameters of the model. 



Parameter 


Symbol 


Value 


Distance 


d 


6 kpc 


Black hole mass 




6M 


Viewing angle 


<t> 


30° 


Jet bulk Lorentz factor 


Tjct 


2 


Jet injection point 


20 


4.5 X 10 7 cm 


Ratio ro/zo 


X 


0.1 


Ratio 2Ljet/iEdd 


'/jet 


> 10" 4 


Ratio L Te \/Lj et 


9rel 


0.1 


Ratio L p /L e 


a 


> 1 


Magnetic field decay index 


in 


1 - 2 


Ratio Us/U^ -m) at z acc 


P 


0.1 - 1 


Particle injection index 


a 


1.5 


Minimum particle energy 




3; 2mc 2 


Acceleration efficiency 


V 


10" 4 - 0.1 



All parameters were kept fixed during fitting except gj et , a, r], 
E mm and p, that varied in the range indicated. 
The value adopted for the jet Lorentz factor is typical for micro- 
quasar jets (e.g. Fender, Belloni & Gallo 2004). 



summarized in Table [2] We allowed Qjot, a, r/, E mm and p 
to vary subject to the constraints discussed above, whereas 
the rest of the parameters were kept fixed during the fitting 
process. We performed least-squares fits (Kay 1993) to the 
observational data. The quality of the fits was quantified 
calculating the correspondent value of x 2 , 

x 2 =J2 (Fo l? Fm)2 ■ ( 45 ) 

obs 

Here F h s is the observed flux, F m is the value predicted 
by the model, and AF D b s is the uncertainty associated with 
every observational point. The best fit for a given set of 
parameters was found minimizing the value of \ 2 ■ 



3.2 Spectral energy distributions 

Figures [2] and [3] show the best fits obtained for a set of 
simultaneous radio and X-ray data taken on February 3rd 
1997 and April 2nd 1999 (Obsl and Obs2 in Tabled] respec- 
tively), when the source was in a luminous low-hard state. 
In the case of Figure [2] the location of the acceleration re- 
gion z&cc was determined through condition (J5J), whereas in 
the case of Figure [3] condition (|8} was applied. 

Each graphic corresponds to a different value of the 
magnetic field decay index m. This parameter strongly de- 
termines the shape of the spectrum, since it fixes the value 
of the field along the jet and consequently z acc . Larger val- 
ues of m yield z aC c closer to the jet base where the magnetic 
field is stronger. The values of the parameters result from 
the fitting are listed in Table [3] 

The spectral energy distributions (SEDs) presented in 
Figure[2]all correspond to Obsl. The best fit is obtained for 
m = 1.2. The X-ray data range is always covered by the 
synchrotron emission of primary electrons, but as m grows, 
synchrotron radiation of secondary pairs begins to dominate 
at radio wavelengths. This diminishes the quality of the fit. 
Also as m increases the slope of the X-ray spectrum gets 
worse modeled, indicating that the injection index should 
be harder than the assumed a — 1.5. 

Between ~ 1 GeV and ~ 1 TeV, the emission is dom- 
inated by synchrotron self-Compton (SSC) radiation and 
synchrotron emission of protons and secondary particles; 
at higher energies the contributions of proton-proton and 
proton-photon interactions are the dominant ones. All these 
processes become more relevant when z acc is nearer the jet 
base, since the magnetic field is stronger and enhances the 
synchrotron radiation of p, tt and fj, . Also the matter 
and photon densities are larger, providing denser targets for 
pp and P7 collisions, and SSC scattering. The contribution 
of secondary pairs from 77 annihilation is significantly in- 
creased for large values of m due to this effect as well. In 
all cases the best fits favour large minimum particle ener- 
gies, E-axin ~ lOOmc 2 . A powerful jet (qj ct ~ 0.8 — 0.9) and 
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Log (E, / eV) Log (E. ( / eV) 

Figure 2. Best fit spectral energy distributions (SEDs) of Obsl for different values of the magnetic field decay index m. The graphics 
are labeled as in Table [3] The position £ acc of the acceleration region was determined demanding that Us < U^. See Tables [2] and [3] for 
the values of the rest of the parameters. The subindices (77), [pj) and (fi) indicate pairs created through photon-photon annihilation, 
photopair production and muon decay, respectively. The thick lines are the sensitivity limits of Fermi and HESS and CTA. 




Log(E r /eV) Log M (E.,/eV) 

Figure 3. The same as in Figure [2] but with £ acc calculated from the condition Ub < U m - Model E is a fit to Obsl and model F to 
Obs2. 
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Table 3. Best fit values of the relevant parameters. 



Model 


Obs 


a 


in 


E min /mc 2 


STjet 


a 


'/ 


p 


■Sacc [-^<?] 


B(z acc ) [G] 




X 2 /d.o.f. 


A 


Obsl 


1.5 


1.2 


97.4 


0.9 


1.5 


0.1 


0.1 


9.7 x 10 4 


8.4 x 10 3 


3.8 x 10 38 


1.42 


B 


Obsl 


1.5 


1.5 


99.2 


0.8 


1.1 


0.08 


0.75 


1.4 x 10 2 


1.5 x 10 7 


8.6 x 10 40 


2.0 


C 


Obsl 


1.5 


1.8 


96.3 


o.s 


1.6 


0.03 


0.5 


1.4 x 10 2 


1.4 x 10 7 


8.4 x 10 40 


3.1 


D 


Obsl 


1.5 


2.0 


92.4 


0.75 


1.4 


0.03 


0.75 


85.5 


2.3 x 10 7 


1.3 x 10 41 


3.3 


E 


Obsl 


1.5 


1.8 


15.0 


1.00 


1.5 


0.1 


0.1 


1 x 10 4 


5.0 x 10 3 


4.5 x 10 40 


0.8 


F 


Obs2 


1.5 


2.0 


11.5 


0.73 


2.7 


0.1 


0.1 


3.5 x 10 3 


1.2 x 10 4 


3.0 x 10 40 


0.98 


G 


Obs3 


1.5 


2.0 


2.0 


6.4 x 10" 3 


2.0 


3 x 10~ 3 


0.4 


3 x 10 3 


3.6 x 10 3 


1.4 x 10 35 


0.15 


H 


Obs4 


1.5 


2.0 


2.0 


6.3 x 10" 3 


2.0 


1 x 10~ 4 


0.15 


4.9 x 10 3 


1.4 x 10 3 


7.5 x 10 34 


0.15 


I 


Obs3 


2.2 


2.0 


25.2 


6.6 x 10~ 3 


2.0 


0.1 


0.4 


2 x 10 2 


3.1 x 10 3 


2.4 x 10 37 


0.8 


J 


Obs5 


1.5 


2.0 


92.7 


1.0 


1.0 


0.1 


0.1 


6 x 10 3 


4 x 10 3 


1 x 10 42 


6.2 


The values of < 


?jot, a, 


■q, E m 


in and p were 


left free during the fitting procedure; 


X 2 /d.o.f. is 


the value of 


X 2 per 





degree of freedom. 



equipartition of energy between primary protons and lep- 
tons (a « 1) is also required, since the power injected in 
electrons needs to be as large as possible to account for the 
X-ray observations. 

Models E and F in Figure ^correspond to fits of Obsl 
and Obs2, respectively. In both cases, z acc was calculated 
demanding that (7b < U m . For the same m and p, this con- 
dition gives larger values of z acc and weaker magnetic fields. 
Now the best fits are obtained for large values of m. These 
sets of parameters reproduce the slope of the X-ray spec- 
trum for the same value of the injection index better than 
the models of Figure [2] 

Figure [4] shows two model fits to low-luminosity low- 
hard state observations of GX 339-4, carried out in 1999 
(Obs3 and Obs4). The radio and X-ray emission is due to 
primary electrons; all radiative contributions of protons and 
secondary particles are negligible. The jet power required to 
account for the data is now only a fraction qj ct « 6 x 10 -3 of 
the Eddington luminosity. The best fit models are obtained 
for low values of the acceleration efficiency and minimum 
particle energy. 

For each model we calculated the synchrotron emission 
of thermal electrons at the base of the jet. For an electron 
energy E e ~ 2m e c 2 and a magnetic field Bo « 10 6-7 G, 
the peak of the spectrum is at _E 7 10 eV. However, the 
luminosity of this component is below or just above the jet 
emission. This contribution is not significant in the relevant 
energy bands. Therefore, the results of the fits are not af- 
fected. 



3.3 Radio/X-rays and NIR/X-rays correlations 

The analysis of simultaneous radio and X-rays observations 
from 1997-1999, led Corbel et al. (2003) to find out that 
the fluxes in both energy bands are tightly correlated. They 
showed that the radio flux at 8.6 MHz is related the 3-9 keV 
integrated X-ray flux as Fr, oc AF^. This correlation sug- 
gests a common origin in the je10 (synchrotron radiation). 
According to Markoff et al. (2003), if all model parameters 



3 Alternative models to explain the radio/X-rays correlation have 
been suggested. Markoff et al. (2005) presented fits to simultane- 
ous radio and X-ray data of GX 339-4 obtained applying a corona 
model. Furthermore, in Heinz & Sunyaev (2003) it is shown that 
for an ADAF-like boundary condition, the radio flux from the 



except the jet power are kept frozen, the correlation index 8 
is given by 



8 = 



17/12 - 2/3<fe 
17/12 - 2/3<5> 



(46) 



Here <5x is the spectral index of the X-ray region of the 



synchrotron spectrum (Fx. 



and 5r that of the syn- 



chrotron radio flux (Fb. oc u Sr ). 

In our model radio and X-ray emission is due to syn- 
chrotron radiation of electrons. We find a value of the radio 
spectral index <5r ~ 0.33, which corresponds to the opti- 
cally thick part of the spectrum from a particle distribution 
with a low-energy cutoff. The value of the X-ray spectral 
index is <5x ~ —0.8, as expected for an injection function 
Q e oc E~ 5 (notice that electrons are then strongly cooled 
due to synchrotron losses). These values yield 8 ~ 0.6. Fig- 
ure [5] shows the correlation curves predicted by our model 
for cases A, E, F and G, together with the data from Corbel 
et al. (2003). The model results are in reasonable agreement 
with the observations. 

Simultaneously with the radio and X-ray observations 
of Obs3 and Obs5, GX 339-4 was also detected at NIR and 
optical wavelengths (Markoff et al. 2003; Homan et al. 2005). 
As in the case of the radio emission, the NIR/optical flux 
also displays a strong correlation with the X-ray flux. 

From an analysis of low-hard state data from the 2002 
outburst of GX 339-4 (same epoch as Obs5), Homan et 
al. (2005) shown that the flux density in the NIR //-band 
and the 3 — 100 keV bolometric X-ray flux correlate as 
Fh oc AFx, with 8 = 0.53. A similar correlation was found 
between the optical V-band and /-band flux densities and 
the integrated X-ray flux, with correlation indices 8 = 0.44 
and 8 — 0.48, respectively. These correlations disappear 
when the source leaves the low-hard state. The i/-band emis- 
sion, however, rises and decays faster than the optical dur- 
ing the state transition, while the slope between the / and 
V bands remains constant. As suggested by Homan et al. 
(2005), this may indicate a different origin for the NIR and 
optical emission during the low-hard state. 

The correlations between the radio/X-ray and NIR/X- 
ray fluxes suggest that the emission in the three ranges must 
originate in the jet. This is further supported by the fact 



base of the jet scales with the black hole mass and the accretion 
rate, independently of the assumed jet model. 
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Log(E y /eV) Log(E y /eV) 

Figure 4. Best-fitting SEDs of Obs3 (model G) and Obs4 (model H). The decay index of the magnetic field is m = 2, and the position 
Zacc of the acceleration region was determined demanding that Ub < Um- See Tables [2] and \3\ for the values of the rest of the parameters. 
The sensitivity limits of Fermi, HESS and CTA are indicated. Optical data in model G (not shown, see Figure [Sj were not included in 
the fit. 



I 1 1 1 I 1 

■ Data 


i 1 i ■ i 


Model A 




■ Model E 




Model F 




- - Model G 




i.i.i. 





-1.5 -1.0 -0.5 0.0 0.5 1.0 

Log 10 AF x (3 - 9 keV) (10 10 erg 1 s 1 cm 2 ) 

Figure 5. Radio/X-ray flux correlations in GX 339-4. The dif- 
ferent curves correspond to models A, E, F and G. In each case, 
the slope was calculated as in Markoff et al. (2003), and to de- 
termine the intercept we used corresponding SEDs in Figures [2] 
Q5] and [3] The model correlation index is S ~ 0.6, whereas that of 
the 1997-1999 observational data is 6 ~ 0.7 (Corbel et al. 2003). 



that the NIR flux extrapolates back to the radio data (see 
also Corbel et al. 2003). Direct or reprocessed emission from 
an accretion disc can be ruled out due to the shape of the 
NIR/optical spectrum and the short decay time scales. Fur- 
thermore, the NIR and radio fluxes are quenched when the 
disc begins to contribute significantly to the X-ray emission. 
Homan et al. (2005) conclude that the NIR emission prob- 
ably originates in the jet, and approximately coincides with 
the position of the break of the synchrotron spectrum. The 
optical flux may be due to thermal and/or non-thermal re- 
processed radiation from the accretion disc or star, or from 
a region of the jet different from where the NIR emission is 
produced. 



These ideas are supported by the recent results of Coriat 
et al. (2009), who presented an analysis of five years of obser- 
vations of GX 339-4 (from 2002 to 2007, a period that com- 
prises five outbursts). They found a strong IR/X-ray correla- 
tion over four decades in flux during the low-hard state. The 
correlation index, however, is not unique: a break appears at 
bolometric (3-100 keV) X-ray fluxes ~ 1.1 x 10" 10 erg s" 1 
cm" 2 (L x ~6x 10" 4 L Ed d for M = 6M© and d = 6 kpc). 
Coriat et al. (2009) argue that this break can be explained 
attributing the X-ray emission to SSC radiation from the 
jet (see also Nowak et al. 2005, where it is suggested that 
models more complex than a single jet synchrotron compo- 
nent maybe needed to explain the correlations) . They find no 
clear evidence of a similar break in the l/-band/X-ray corre- 
lation, and suggest that the optical emission in the low-hard 
state is dominated by the outer part of the accretion disc, 
and not by the jet. 

This correlation is not peculiar of GX 339-4, but it 
seems to be a signature of low-mass black hole X-ray bina- 
ries. Russell et al. (2006) analysed radio, NIR, optical and 
X-ray data from 16 sources (including extragalactic systems 
in the LMC). Their results agree with those of Homan et al. 
(2005) for GX 339-4. They estimate that the jet contribution 
to the NIR emission during the low-hard state is ~ 90%, but 
only ~ 50% to the I and V bands. 

Figure [S] shows the best fit models obtained for Obs3 
(including the lowest energy point in the optical) and Obs5. 
In the case of Obs5, the radio, NIR, optical and X-ray data 
are reasonably well reproduced using a hard particle injec- 
tion spectral index a — 1.5. In the case of Obs3, the data 
at optical frequencies cannot be accounted for with a single 
synchrotron component. However, using a softer particle in- 
jection spectral index a — 2.2, it is possible to obtain models 
models where the synchrotron turnover occurs in the opti- 



Leptonic/hadronic models for MQs 11 




Log (Ey / eV) Log (E^ / eV) 



Figure 6. Best fit SEDs of Obs3 (model I) and Obs5 (model J). In model I, the arrow indicates the only point in the optical band that 
was included in the fit. The decay index of the magnetic field is m = 2. The position z aC c of the acceleration region was determined 
demanding that Ub < in Model I, and Ub < U m in Model J. See Tables [2] and [3] for the values of the rest of the parameters. The 
sensitivity limits of Fermi, HESS and CTA are indicated. 
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Figure 7. Attenuation factor for a photon of energy due to 
77 annihilation in the total photon field of the emission region. 
Absorption is negligible in all cases and does not modify the pro- 
duction spectrum. 

caQ. The rise in the spectrum at higher energies, however, 
cannot be fitted. This emission must have a different origin, 
for example in an accretion disc (Markoff et al. 2003). 

3.4 Photon absorption 

In order to asses the effect of photon self-absorption by 
77 annihilation, we calculated the attenuation parameter 
exp(— t 77 ) in equation (|43[l . As it can be seen from Figure 
[7] internal attenuation is almost negligible due to the low 
photon density in the emission region. The production spec- 
trum is not appreciably modified in any case, contrary to 
some models calculated by Romero & Vila (2008). 

4 Synchrotron radiation of thermal electrons from the base of the 
jet is not relevant in these models either. 



3.5 Discussion 

According to the results of Figures [5] EH and [5] during the 
high luminosity low-hard state, GX 339-4 would be a source 
detectable by Fermi in the energy range 100 MeV - 10 GeV. 
In this range the emission is due to electron synchrotron ra- 
diation and SSC. However, it will be far from being a bright 
source, possibly appearing with a significance of 5-6 a af- 
ter one year of exposure. At higher energies, the emission 
of hadronic origin will be harder to detect with the present 
instruments. For some models, the predicted luminosities at 
~ 1 TeV are above the sensitivity of HESS for a point source 
at 6 kpc through exposures of more than 50 hours. Future 
instruments like the Cherenkov Telescope Array (CTA) will 
make it easier to obtain a firm detection. Such a detection at 
high energies is crucial to evaluate the hadronic content of 
the jet. On the other hand, when the source is in a low lumi- 
nosity phase, according to the results in Figure [Jit would be 
undetectable at high energies with the present instruments. 



4 POSITRON PRODUCTION RATE 

Recently, an analysis of several years of data collected by the 
INTEGRAL satellite led Weidenspointner et al. (2008) to 
suggest that the observed asymmetry in the distribution of 
the electron-positron annihilation line traces the distribution 
of hard low-mass X-ray binaries (LMXBs) in the Galaxy. 
These authors discovered that the 511 keV line emission 
from the inner galactic disc is clearly asymmetric: the flux 
from the region of negative galactic longitudes is 1.8 times 
larger than that from the corresponding region of positive 
longitudes. The same type of asymmetry is observed in the 
spatial distribution of hard LMXBs, those that show ap- 
preciable emission above 20 keV. The authors argue that 
positrons escaping the injection region with an energy of 
~ 1 MeV will not diffuse more than ~ 100 pc before an- 
nihilating. This distance is short enough for the line emis- 
sion to be still correlated with the large-scale distribution 



12 G. S. Vila and G. E. Romero 



of the sources. Weidenspointner et al. (2008) estimate that 
the positron production rate required to account for the ob- 
served flux is ~ 10 41 erg s _1 . 

To investigate this possibility we calculated the positron 
injection rate in our model. According to Heinz (2008), the 
number of injected positrons per unit time can be roughly 
estimated as 



N e+ 



(47) 



In this expression L e ± is the total luminosity injected in 
pairs and 7 is the mean Lorentz factor of the positrons when 
they leave the source. 

It is reasonable to expect that positrons have almost 
completely cooled when they reach the end of the jet, 
and thus 7 is of the order of the jet bulk Lorentz factor, 
7 ~ Tjct = 2. In our case, the most relevant process of pair 
production is 77 annihilation. The calculated positron in- 
jection rates are shown in Table For those models that 
correspond to the high X-ray luminosity states, N e + is very 
near or even exceeds the necessary lower limit given by Wei- 
denspointner et al. (2008). There are about > 100 LMXB in 
the Galaxy and, although not detected yet, possibly most of 
them produce jets. Even if many of them are considerably 
less powerful than the jet in GX 339-4, the added contribu- 
tions might well account for the observed flux at 511 keV. 
Our estimations show that the proposed association between 
LMXB and the annihilation line emission is indeed feasible 
at least in energetic terms. In this way, there may be no 
need to resort to other more exotic explanations, such as 
annihilation of dark matter. A more detailed discussion on 
this topic will be presented elsewhere. 



5 SUMMARY AND CONCLUSIONS 

We have shown that, under certain general conditions, the 
model developed here is capable of explaining the observed 
radio and X-ray spectrum of the low-mass microquasar GX 
339-4. The parameter that mainly determines the character- 
istics of the obtained spectra is the magnetic field decay in- 
dex m. This parameter, in turn, depends on unknown factors 
as the geometry and turbulence level of the magnetic field 
in the inner jet. The best fits to the high luminosity 1997 
and 1999 data are obtained in those cases where the acceler- 
ation/emission region is placed relatively far from the com- 
pact object and the jet apex, z aC c ~ 10 9 and z acc ~ 10 11 cm 
(models A and E, respectively). In both cases the value of the 
magnetic field is similar, B(z acc ) ~ 10 4 G, although m = 1.2 
in model A and m = 2 in model E. The slope of the X-ray 
spectrum was difficult to reproduce, even assuming a hard 
particle injection with spectral index a = 1.5. Nevertheless, 
a harder particle injection could be achieved in principle 
through diffusive acceleration at relativistic shocks, see for 
example the works of Stecker, Baring & Summerling (2007) 
and Keshet & Waxman (2005). In any case, the high X-ray 
luminosities require a powerful jet with a large leptonic con- 
tent. In fact, all the fits yield a ~ 1, what means that as 
much energy is given to the primary relativistic electrons as 
it is allowed by the constraints imposed. The hadronic con- 
tribution to the spectrum in cases A and E is undetectable 
with the present gamma-ray instruments. In the other mod- 
els, synchrotron radiation of protons and secondary muons 



and pions, and at higher energies the contribution of pp in- 
teractions, could be detectable by Fermi and HESS (and in 
the future by CTA), respectively. 

For the low-luminosity observations (models G and H), 
the best fits were obtained for low values of the acceleration 
efficiency and minimum particle energy. The required jet 
luminosity is now only Lj ct ~ 6 x 10 -3 I/Edd- The predicted 
emission above ~ 100 MeV is too faint to be detected with 
the present gamma-ray telescopes. 

We have also calculated fits to simultaneous radio, 
NIR/optical and X-ray observations from 1999 and 2002 
(models I and J). For these sets of parameters, the break 
in the synchrotron spectrum occurs approximately in the 
NIR, and the lowest-energy data was reasonably fit. The 
rising shape of the spectrum at optical wavelengths, how- 
ever, could not be reproduced. This component is likely to 
originate mostly outside the jet, probably in the accretion 
disc. 

In all models, the spectrum is essentially of leptonic 
origin. In this sense, the results do not differ from those of 
previous works like those of Markoff et al. (2003, 2005). Our 
model, however, besides making predictions for the emission 
in the very high-energy regime, introduces some refinements 
over the previous scenarios adopted for this source. The par- 
ticle distributions are calculated self-consistently taken into 
account the effect of energy losses on the injection spec- 
trum. We also calculate the radiation emitted by secondary 
particles produced in hadronic interactions, and that of the 
electron-positron pairs product of photon-photon annihila- 
tion. The importance of photon self-absorption is assessed 
as well, although it turns out not to be relevant since the 
emission region is in a zone of low radiation density. 

We have also shown that the pair injection rate is signif- 
icant enough, if this kind of model is solid in general for low- 
mass microquasars, to account for the observed line emission 
at 511 keV, according to the lower limit given by Weidens- 
pointner et al. (2008). If the proposed association between 
hard low-mass X-ray binaries and the electron-positron an- 
nihilation line flux can be proved, other explanations such 
as annihilation of dark matter could result unnecessary. De- 
tailed models for dealing with this topic are in preparation. 
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